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Abstract 

The generalization of dynamical scaling to local scale invariance is 
reviewed. Starting from a recapitulation of the phenomenology of 
ageing phenomena, the generalization of dynamical scaling to local 
scale transformation for any given dynamical exponent z is described 
and the two distinct types of local scale invariance are presented. The 
special case z = 2 and the associated Ward identity of Schrodinger 
invariance is treated. Local scale invariance predicts the form of the 
two-point functions. Existing confirmations of these predictions for (I) 
the Lifshitz points in spin systems with competing interactions such 
as the ANNNI model and (II) non-equilibrium ageing phenomena as 
occur in the kinetic Ising model with Glauber dynamics are described. 
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1 Ageing and dynamical scaling 



Non-equilibrium critical phenomena are a subject of intense research activity 
in physics. Rather than providing a general and exhaustive definition, we 
shall present a typical example of current interest, chosen such that our 
main question can be introduced in a natural way. 

A common way to reach a situation of non-equilibrium criticality is through 
a rapid change of one of the macroscopic variables which enter into the equa- 
tion of state. For definiteness, consider a simple ferromagnet. The simplest 
model system of this kind is the Ising model which considers magnetic mo- 
ments attached to the sites of a (regular) lattice A. By hypothesis, these 
moments may only take two values, say 'up' and 'down', and are described 
through a spin variable ai = ±1 attached to the lattice site i G A. In equi- 
librium, to a given configuration of spins {a} = {cTi}i^\ one associates an 
energy 7-^{cr} = — ^i^j where the sum is usually restricted to nearest- 

neighbour pairs on the lattice. The probability of a given configuration {cr} 
in equilibrium is then given by 



where T is the temperature. It is well-known that systems of this kind show 
an order- disorder phase transition, at least m. d > 1 dimensions, such that 
the order parameter := (^i^^cri) has a non- vanishing value for tempera- 
tures below a certain critical temperature > and vanishes for T > T^. 
Such a model may be brought out of equilibrium, for example, by starting 
initially from a fully disordered state (effectively corresponding to an infinite 
temperature) and then quench the system rapidly to a fixed temperature 
below the system's critical temperature Tc > 0. Empirically, it is then found 
that the system evolves slowly towards the equilibrium state corresponding 
to the fixed temperature T. However, the details of this relaxation process 
may be complex. In many materials undergoing such thermal treatment, the 
time-dependent properties may well depend on the entire thermal history of 
the sample. The resulting ageing behaviour has been in the focus of intensive 
study, see [H 12 13 HI E] for reviews. 

In the Ising model example, a popular choice to simulate the update of 
the individual spins is through Glauber dynamics [Hj. Glauber dynamics may 
be realized through the discrete-time heat-bath rule ai{t) ^ ai{t + 1) such 
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Figure 1: Snapshots of the coarsening of ordered domains in the 2D Glauber- 
Ising model, after a quench to T = 1.5 < from a totally disordered state 
and at times t = 25 (left) and t = 275 (right) after the quench. 

that 

ai{t + 1) = ±1 with probability | [1 ± tanh(/ii(t)/T)] (2) 

with the local field hi{t) = where runs over the nearest neigh- 

bours of the site i. The main property of this dynamical rule is that the order 
parameter is not conserved. For non-equilibrium systems, there is no gen- 
eral theoretical framework for the calculation of the time-dependent probabil- 
ities P({(j}; t). However, with the choice P({cr}; t) can be found exactly 
in ID by solving an associated master equation [H]. The time- dependent 
spin-spin correlators and their approach towards equilibrium can be deter- 
mined exactly in this special case. Otherwise, numerical simulation must be 
used. 

In figure ^ W6 illustrate the evolution of the microscopic state after a 
quench to T < T^. We see that ordered domains form and slowly grow, with 
a typical size L{t) ~ t^^^, where z is the dynamical exponent. The presence of 
algebraic growth laws for such dynamic quantities, although the equilibrium 
state of the model is not critical, is characteristic for the ageing phenomena 
we wish to consider. Indeed, inside a given domain the spins are completely 
ordered and the time evolution of the model only occurs through the slow 
motion of the domain walls. 

This evolution is more fully revealed through the study of two-time quan- 
tities, such as the two-time correlation function C(t,s;r) and the (linear) 
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Table 1: Non-equilibrium exponents a and b for quenches from a fully disor- 
dered state onto and below the critical point Tc > of simple ferromagnets, 
according to [7j. 
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h=0 



where is the order parameter, h the conjugate magnetic field, t is called 
the observation time and s the waiting time. Space translation invariance 
is assumed throughout. Autocorrelators are given by C{t,s) = C(t, s;0) 
and autoresponses by R{t, s) = R{t, s; 0).^ One says that a system is ageing 
when C{t,s) and/or R{t,s) depend on both t and s and not only on their 
difference. Ageing phenomena are therefore related to a breaking of time- 
translation invariance. Physically, ageing occurs in the regime when s and 
T = t — s > are simultaneously much larger than any microscopic time scale 
Tmicro- In many systems, one finds in the ageing regime a dynamical scaling 
behaviour, see jSl Ej 

s) = s~'fcit/s) , Rit, s) = s-'-'^fRit/s) (4) 

where a and b are non-equilibrium exponents. We emphasize that this non- 
equilibrium scaling also occurs for quenches exactly to the critical point 
T = Tc, but the values of the exponents found at T = will be different 
from those found for T < Tc. For T < Tc, it can be shown from phenomeno- 
logical arguments that z = 2 for dynamical rules with a non-conserved order 
parameter while at T = Tc, z has a non-trivial value which must be 
determined from a full renormalizat ion-group approach. 



■^We implicitly assume here that the mean order parameter — 0, otherwise C{t, s) 

must be replaced by the connected correlator r{t, s) = C{t, s) — {(f>{s)). 
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In table Q we collect the values of the non-equilibrium exponents a and 
h. Indeed, it turns out that for T < the equilibrium behaviour of the 
order parameter correlator is essential. If Ceq(r) ~ e~'^'/^ with a finite ^, we 
say that the system is of class S, while if Ceq(r) ~ |r-|~('^~2+''), the system 
is said to be of class L, where r/ is a standard equilibrium critical exponent. 
The results for a are based on the well-accepted physical image that ageing 
comes from the slow motion of the domain walls [21 El E]- This idea had 
been questioned recently and it was argued that because of some anomalous 
scaling behaviour, the simple relation a = 1/2 should be invalid in the 2D 
Glauber-Ising model (which belongs to class S) and rather a = 1/4 here pUj . 
However, recent tests have refuted this claim and reestablish that indeed 
a = 1/2 in this model E]. 

The scaling functions are expected to behave for large arguments x = 
t/s^ 1 asymptotically as 

fc{x) ~ x""^/^ , fR{x) ~ x-^^'^ (5) 

where Ac and \r are called the autocorrelation exponent [121112] and autore- 
sponse exponent ^3], respectively. For a fully disordered initial state, it is 
traditionally accepted that Xc = Xr = X. However, for spatially long-ranged 
initial correlations of the form Cini(r) ~ (with a < 0) the relation 

Xc = Xr + u (6) 

has been conjectured ^3] and Ac = Xr is only recovered if a renormalizes 
to zero.^ Furthermore, the rigorous arguments of [18j readily yield A^ > 
{d + a)/2, for simple ferromagnets without disorder. Very recently, distinct 
exponents Ac 7^ Xr have also been found in the random sine-Gordon model 
and in addition A^ < d/2 violates the rigorous bound mentioned above (TU]. 

Another central question in this context is whether/when under the con- 
ditions just described the system is in thermodynamic equilibrium. It is 
convenient to consider the fluctuation-dissipation ratio [201 12I| 

X{t,s)=TR{t,s)(^^^^^ ' (7) 

^The conjecture ^ is in agreement with all known results from the kinetic spherical 
model |14lll5j and also explains those of the 2D XY model with T < Tkt jIQ- Remarkably, 
for the Ising model with long-range initial conditions, only Ac appears to be known |17j . 
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At equilibrium, the fluctuation-dissipation theorem states that X{t,s) = 1. 
The breaking of the fluctuation-dissipation theorem has been investigated 
intensively both theoretically (see e.g. [3111122112301211) and experimentally 

[23l2nil2Il. 

Having thus reviewed the main phenomenological aspects of the ageing 
behaviour of simple ferromagnets, we can now formulate our main question: 
is it possible to generalize the dynamical scaling described by eq. ^ to a 
local scale invariance ? By this we mean that we wish to generalize the 
transformations of dynamical scaling t ^ Vt and r i— > 6r to a local form 
where h = b{t,r) |2H1- We therefore seek to extend dynamical scaling in 
a way analogous to the extension of global scale invariance in equilibrium 
critical phenomena to conformal invariance. 

This review is organized as follows. In section 2 we present the main 
points of the construction of local scale transformations, for an arbitrary 
dynamical exponent z. In section 3 we consider in more detail the case 
z = 2 where the relation to the conformal group and local Ward identities 
are discussed. In section 4, general expressions for the scaling functions of 
two-point functions are derived. Application to Lifshitz points are briefly 
discussed in section 5. Finally, in section 6 we return to applications of local 
scale invariance to ageing systems, focussing in particular on the prediction 
of the scaling of the response functions and tests in the Glauber-Ising model. 

2 Construction of local scale transformations 

We are interested in systems with strongly anisotropic or dynamical criti- 
cality. By definition, two-point functions of such systems satisfy the scaling 
form 

G{t, r) = h^^'GihH, hr) = [rt-^l^] = r^^^^n (tr^^) (8) 

where t stands for 'temporal' and r for 'spatial' coordinates, X IS db scaling 
dimension, 6 the anisotropy exponent (when t corresponds to physical time, 
6 = z is called the dynamical exponent) and Q are scaling functions. 
Physical realizations of this are numerous, see jT| |21 El El 1211 and references 
therein. For isotropic critical systems, 6 = 1 and the 'temporal' variable 
t becomes just another coordinate. It is well-known that in this case, scale 
invariance ^ with a constant rescaling factor b can be replaced by the larger 
group of conformal transformations b = b{t, r) such that angles are preserved. 
It turns out that in the case of one space and one time dimensions, conformal 
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invariance becomes an important dynamical symmetry from which many 
physically relevant conclusions can be drawn |3(Jj . 

Given the remarkable success of conformal invariance descriptions of equi- 
librium phase transitions, see e.g. [^|32j for introductions, one may wonder 
whether similar extensions of scale invariance also exist when Q Indeed, 
for 6 = 2 the analogue of the conformal group is known to be the Schrodinger 
group [SniElj (and apparently already known to Lie). We shall first describe 
the construction of these local scale transformations for arbitrary 9 ^ 1,2, 
show that they act as a dynamical symmetry, then derive the functions $, Vt 
and finally comment upon some physical applications. We shall present the 
main results as formal propositions and refer to |2H| for details and the proofs. 

The defining axioms of our notion of local scale invariance from which 
our results will be derived are as follows (for simplicity, va. d = 1 space 
dimensions) 



1. We seek space-time transformations with infinitesimal generators 
such that time undergoes a Mobius transformation 

t^t' = ^^^ a5-/37 = l (9) 
-ft + 6 ^ ' ^ ^ 

and we require that even after the action on the space coordinates is 
included, the commutation relations = {n — m)Xn+m remain 

valid. This is motivated from the fact that this condition is satisfied 
for both conformal and Schrodinger invariance. 

2. The generator Xq of scale transformations is Xq = —tdt — 9~^rdj. —x/9 
with a scaling dimension x. Similarly, the generator of time translations 
is X_i = -dt. 

3. Spatial translation invariance is required. 

4. Since the Schrodinger group acts on wave functions through a projective 
representation, generalizations thereof should be expected to occur in 
the general case. Such extra terms will be called mass terms. Similarly, 
extra terms coming from the scaling dimensions should be present. 

5. The generators when applied to a two-point function should yield a 
finite number of independent conditions, i.e. of the form X„G = 0. 
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Proposition 1: Consider the generators 



X„ = -t-^% - E (I t 1 A,or'''-H-'^dr - E U + 1 ^'^o^''^'^"' 
fc=0 ^ ^ fc=o ^ ^ 

(10) 

where the coefficients A^o and B^o are given by the recurrences An+i.o = 
OAnoAio, Bn+ifi = (ra5„ov4io - A„ofiio) for n > 2 where Aqq = 1/9, 
Bqo = x/6 and in addition one of the following conditions holds: (a) A20 = 
OAlo (h) Aio = A20 = (c) A20 = = (d) Aio = 5io = 0. These 
are the most general linear first-order operators in dt and dr consistent with 
the above axioms 1. and 2. and which satisfy the commutation relations 
[X„, Xni] = [n — n')Xn+n' for all n, n' G Z. 

Closed but lengthy expressions of the X„ for all n G Z are known [28j. In 
order to include space translations, we set ^ = and use the short-hand 
Xn = —t^^^dt — ttndr — bn- We then define 

where m = — y + A; and k is an integer. Clearly, F_Ar/2 = —dr generates 
space translations. 

Proposition 2: The generators Xn and Ym defined in eqs. illOfll\) satisfy 
the commutation relations 

[Xn, = (n- n!)Xn+n' , [Xn, Ym] = (^n^ - rnj Yn+m (12) 

in one of the following three cases: (i) Biq arbitrary, Aiq = A20 = -B20 = 
and N arbitrary, (ii) Biq and B20 arbitrary, Aiq = A20 = and N = 1. (Hi) 
Aio and Biq arbitrary, A20 = AIq, B20 = |Aio-Bio and N = 2. 
In each case, the generators depend on two free parameters. The physical 
interpretation of the free constants Aio, A20, -Bio? -B20 is still open. In the 
cases (ii) and (iii), the generators close into a Lie algebra, see for details. 
For case (i), a closed Lie algebra exists if Biq = 0. 

Turning to the mass terms, we now restrict to the projective transforma- 
tions in time, because we shall only need those in the applications later. It 
is enough to give merely the 'special' generator Xi which reads for Bio = 
as follows [2H1 

X, = -t^dt-Ntrdr-Nxt-ar^dt'-Pr^dl^"-'^"' -^dl^^-'^'^'r^ (13) 



7 



where a,^,^ are free parameters (the cases (ii,iii) of Prop. 2 do not give 
anything new). Furthermore, it turns out that the relation [Xi, Yi\f/2] = for 
N integer is only satisfied in one of the two cases (I) /3 = 7 = which we call 
Type I and (II) a = which we call Type II. The distinction between Type 
I and Type II is very important for the physical applications. 

In both cases, all generators can be obtained by repeated commutators 
of X_i = —dt, Y_j^/2 = —dr and Xi, using (fT^ . Commutators between two 
generators Ym are non-trivial and in general only close on certain 'physical' 
states. One might call such a structure a weak Lie algebra. 

Technically, these results depend on the construction of commuting fractional 
derivatives satisfying the rules 



The standard Riemann-Liouville fractional derivative is not commutative, see e.g. 
j36j . But if one uses the Hadamard definition of the Riemann-Liouville fractional 
derivatives, it is possible to generalize it by adding terms containing derivatives of 
the Dirac delta function. Then the relations ()14() can be proven (28) . 

For N = 1, the generators of both Type I and Type II reduce to those 
of the Schrodinger group. For N = 2, Type I reproduces the well-known 
generators of 2D conformal invariance (without central charge) and Type II 
gives another infinite-dimensional group whose Lie algebra is isomorphic to 
the one of 2D conformal invariance [28j. The physical meaning of this new 
realization remains to be understood. 

Dynamical symmetries can now be discussed as follows, by calculating 
the commutator of the generalized 'Schrodinger-operator' S with Xi. We 
take d = 1 and Bio = for simplicity (it is trival that S commutes with the 
generators X_i = —dt and Y_n/2 = ~dr of time and space translations). 

Proposition 3: The realization of Type I sends any solution ip{t,r) with 
scaling dimension x = 1/2 — {N — 1)/N of the differential equation 



into another solution of the same equation. 

Proposition 4: The realization of Type II sends any solution ip{t.,r) with 
scaling dimension x = {9 — 1) /2 + {2 — 9)'y / {(3 + •y) of the differential equation 



\a—l 



(14) 




(15) 




(16) 
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into another solution of the same equation. 

In both cases, iS is a Casimir operator of the 'Gahlei'-subalgebra generated 
from X-i,Y^N/2 and the generahzed Gahlei-transformation Y-N/2+1. The 
equations ()15I16|) can be seen as equations of motion of certain free field 
theories, where x is the scahng dimension of that free field ip. These free 
field theories are non-local, unless N or 6 are integers, respectively. Of course, 
the known dynamical symmetries of the free Schrodinger and Klein- Gordon 
equations are recovered for = 1 and N = 2. 

Phenomenological consequences will be discussed in section 4. 



3 The case z = 2 and the Schrodinger group 

The origin of local scale invariance may be understood in more detail in the 
special case z = 2. We consider the diffusion equation, in d space dimensions 

{2Mdt - dr ■ dr) (pit, r) = (17) 

For fixed Ai, the Schrodinger group is the maximal kinematic invariance 
group [33j on the space of solutions of eq. (fT7|) . It is defined by the space- 
time transformations [TZ is a rotation matrix and a, (3,^,6 E M) 

, , at + (3 , TZr + vt + a 

1 1 — ^t = , r I — > r = — — ; a5 - f3-f = I 18 

7^-1-0 7^-1-0 

and acts projectively on the solutions 4>it,r) We now present several 

new results on the diffusion equation and the Schrodinger group, again in 
the form of propositions. For details and the proofs we refer to |37j . 

Let sci)^ be the Lie algebra of (fTHjl and conf^ is the complexified Lie 
algebra of the conformal group in d dimensions. Time-translations occur 
in scf)^ and are parameterized by f3. We now treat the 'mass' A4 not as a 
constant but as another variable, see |2H]- The wave function 4>{t,r) which 
solves eq. (fT7|) is replaced hy ip = ipi^Ctyr) through 

0(t,r) = ^ /dCe-'^^^(C,t,r) (19) 



which in turn satisfies the equation of motion 



o2 q2 

2i7rrT^ + 7rT)^(C,i,^)=0 (20) 



9 



While sci)^ acts projectively on the wave functions r) (the phase factor 
depends on A4), the action on ip{(,t,r) is via a true representation. Fur- 
thermore, we have 

Proposition 5: The Lie algebra of the maximal kinetic invariance group of 
equation is On the space of wave functions ip defined through 

eq. one has the embedding 



It is now of interest to study the subalgebras of conf^_,_2- particular, the 
parabolic subalgebras of conf^_,_2 can be classified |37, and we obtain several 
new subalgebras, called age or alt. For the case of one spatial dimension 
c? = 1, we illustrate in figure El their definition through the root space of 
confg = B2. To each point of the diagram corresponds a generator of confg. 
For example, the generator of dilatations Xq is one of the two generators 
sitting at the origin (the Cartan subalgebra f)), the generator X_i of time 
translations has the coordinates (—1, —1) and sits in the lower left corner of 
figure Et'bc and the generator Y_i/2 of space translations has the coordinates 
(0, —1). Subalgebras of confg are identified by considering all convex subset 
in this diagram. 

Two of these parabolic subalgebras (figure |2l3c) still contain the generator 
for the dilatations t —>■ b'^t, r br (which is in the Cartan subalgebra of 
confg) but do not contain time-translations anymore. In particular, age^ and 
alti are candidates for a dynamic symmetry algebra of ageing systems. 

While we have worked here in units such that the 'speed of light' c = 1, 
it is possible to reintroduce c exphcitly into the equations of motion such as 
eq. (fT7|). It had been claimed in the past that sci)^ could be obtained from 
a group contraction from conf,^.,.]^ in the 'non-relativistic limit' c ^ 00 jH^. 
This limit procedure can be made precise in terms of the wave function ()19p 
which satisfies a. {d + 2)-dimensional massless Klein-Gordon equation ()2U|). 
We obtain the result 

Proposition 6: The non-relativistic limit c — *• cxo 0/ the conformal group 
acting on the wave function ip defined through M9(l leads to the map 



Although afii and 5d]i have the same dimension, they are not isomorphic, 
see figure l^fc^c. 



scf)rf C confrf+2 



(21) 



confg olti 



(22) 
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Figure 2: Root space of the complexified conformal Lie algebra confg, in- 
dicated by the full and the open points. The double circle in the center 
denotes the Cartan subalgebra. The generators belonging to the three non- 
isomorphic parabolic subalgebras [HZI are indicated by the full points, namely 
(a) scf)i, (b) dgti and (c) olti. 

The dynamical symmetries of the free diffusion equation may be related 
to the symmetries of the action S of the associated free-field theory. This ac- 
tion may be constructed in the standard fashion following the lines of Martin- 
Siggia-Rose theory, see e.g. [23 and references therein. Working with the 
wave function = (j){t,r), we introduce a compound vector e = {5t,5r) for 
small coordinate changes. Let rj be the phase which arises in the transfor- 
mation of the wave function under local coordinate transformations. For 
local theories we expect the following transformation of the action 



where T^^ is the conserved energy-momentum tensor and is the conserved 
probability current. The second term contains the contributions from the 
initial conditions. This local form permits the derivation of Ward identities 
and we obtain 

Proposition 7: Consider an ageing system described by a local action S 
transforming according to \2'J\) and which is invariant under spatial transla- 
tions, phase shifts, Galilei transformations and dilatations with z = 2. Then 
55* = under the action of the special Schrddinger transformation Xi.^ 
This illustrates the conceptual importance of the requirement of Galilei in- 
variance, besides dynamical scaling with 2 = 2, for full local scale invariance 

^Xi is the generator of the 'special' transformation parameterized by 7 in eq. (|18|l . 




(23) 
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to hold. We stress that time-translation invariance is not required. 

Finally, we mention that there exist infinite-dimensional Lie algebras 
which contain scf)^ as subalgebras. For example, the Schrodinger group (fTHjl 
is a subgroup of the group defined by the transformations t —>■ t' and r ^ r' 
where jHZI 

t' = (3{t) , r' = r\l (3{t) or else t' = t , r' = r - cx{t) (24) 

and j3 and ol are arbitrary functions. Whether this has a bearing on the 
non-equilibrium behaviour of spin systems is still open. 



4 Scaling functions 

From a physical point of view, the wave equations ()15I16|) discussed in section 
2 suggest that the applications of Types I and II are very different. Indeed, 
eq. (fT3j) is typical for equilibrium systems with a scaling anisotropy intro- 
duced through competing uniaxial interactions. Paradigmatic cases of this 
are so-called Lifshitz points which occur for example in magnetic systems 
when an ordered ferromagnetic, a disordered paramagnetic and an incom- 
mensurate phase meet (see [lU] for a recent review). On the other hand, 
eq. (fT^ is reminiscent of a Langevin equation which may describe the tem- 
poral evolution of a physical system, with a dynamical exponent z = 9. In 
any case, causality requirements can only be met by an evolution equation 
of first order in dt- 

We now give the scaling functions $, f2 in eq. (jS| from the assumption 
that the two-point function G transforms covariantly under local scale trans- 
formations. 

Proposition 8: [28j Local scale invariance implies that for Type I with 
6 = 2/N , the function VL{v) must satisfy 

{ad^-^ - v^d, - Nx) ^l{v) = (25) 

together with the boundary conditions Q{0) = Qq and fl{v) ~ flooV~'^^ for 
V ^ oo. For Type II, we have 

{du + e{l3 + 7)wC + 2^(2 - ^)7d ^u) = (26) 

with the boundary conditions $(0) = $o one? $(m) ~ $00"^"^^' for u ^ 00. 
Here ^0,00 and $0,00 are constants. The ratio /5/7 turns out to be universal 



and related to x. From the linear differential equations ()25l26p the scaling 
functions VL{y) and $(m) can be found explicitly using standard methods |2H]- 
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5 Application to Lifshitz points 



Given these explicit results, the idea of local scale invariance can be tested 
in specific models. Indeed, the predictions for Q{v) coming from Type I with 
= 4 nicely agree with cluster Monte Carlo data for the spin-spin and 
energy-energy correlators of the 3D ANNNI model at its Lifshitz point. The 
ANNNI model is defined in terms of Ising spins o" = ±1 on a hypercubic 
lattice. The Hamiltonian is 

H = — ^ a^yz \(^xy{z+l) + (^x{y+l)z + <^{x+l)yz\ + ^xyz ^xy{z+2) 

(27) 

xyz xyz 

where k is a constant. Extensive simulations using a generalization ^TJ 
of the Wolff algorithm |43j and a new method ^3] to reduce finite-size effects 
in the estimates of correlation functions led to a considerable improvement 
in the precision of the estimation of the location of the Lifshitz point, viz. 
kl = 0.270(4) and Tl = 3.7475(50) gU. The Lifshitz-point exponents ol = 
0.18(2), /5l = 0.238(5), 7l = 1.36(3), respectively, of the specific heat, the 
magnetization and the susceptibility were obtained [H].^ These values are 
also in good agreement with the results of a careful two-loop study of the 3D 
ANNNI model within renormalized field theory which gives = 0.160, jSi = 
0.220 and 7l = 1.399 In figure we show the scaling function \E'(t>) : = 
v''fl{v) of the spin-spin correlator. Here ( = {2/9){d± + 9)/{2 + 'Jl/Pl), 
where ci_L = 2 is the number of transverse dimensions without competing 
terms and 6 1/2 within the numerical accuracy of the data (see |^ I45| OS] 
for a fuller discussion of this point). We see the expected collapse of the 
data which establishes scaling. In the inset, we compare the data with the 
prediction of local scale invariance. U 6 = 2/N = 1/2, we have A^ = 4 and 
expect from proposition 8 that Q{v) satisfies 

a - v^—^ - 4xvn(v) = 28) 

dv^ dv 

The are two independent solutions which satisfy the required boundary con- 
ditions, such that the form of Q{v) depends on a single universal parameter 
p. Its value can in turn be estimated by considering universal ratios of mo- 
ments of f2(f), see jH] for details. In the inset of figure El we find perfect 
agreement with the prediction of local scale invariance. 

^The fact that the scahng relation + 2(3^ + 7l = 2 is satisfied up to « 0.8% gives 
an a posteriori estimate of the quahty of the numerical data. 
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Figure 3: Scaling function of the spin-spin correlator at the Lifshitz 

point in the 3D ANNNI model. In the main plot the scaling of selected data 
for several values of r± is shown. In the inset, the full set of numerical data 
(gray points) is compared to the solution of eq. (^Hj). The data are from jUj. 



Replacing the Ising model spins by so-called spherical spins G M such 
that YlreA — where M is the number of sites of the lattice, one obtains 
the exactly solvable ANNNS model. At its Lifshitz point, the spin-spin corre- 
lator agrees exactly with the prediction of local scale invariance [lUJ 123123 

On the other hand, the predictions of Type II have been tested extensively 
in the context of ageing ferromagnetic spin systems which are going to be 
described in the next section. 



6 Application to ageing 

Time translation invariance does not hold in ageing systems. The simplest 
way to do this is to remark that the Type Il-subalgebra spanned by Xq^Xi 

^The ANNNS model may be generalized such that Lifshitz points of second order 
appear as the endpoints of lines of Lifshitz points. This allows for a test of local scale 
invariance for = 6 j28| . 
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and the leaves the initial line t = invariant, see (jl3|) . Therefore the 
autoresponse function R(t, s) is fixed by the two covariance conditions XqR = 
XiR = 0. Solving these differential equations leads to 

Proposition 9: jl?! |2H] For a statistical n on- equilibrium system which sat- 
isfies local scale invariance the autoresponse function R{t, s) takes the form 

i?(t,s) = ro(t/s)'+'^-^^/"(t-s)"'"'^ , t>s (29) 

where a and Xr/ z are the non- equilibrium exponents defined in section 1 and 
To is a normalization constant. 

The functional form of R is completely fixed once the exponents a and Xr/z 
are known. Similarly, R{t,s]r) = s)$ (r(t — gives the spatio- 

temporal response, with the scaling function determined by (j26|) j2H]- 
In the z = 2 special case the above result takes a particularly simple form. 
Proposition 10: |48| l^ For an ageing system with dynamical exponent z = 
2 and which is invariant under agc^, the two-time spatio-temporal response 
function is 



R{t, s; r) = R{t, s) exp 



M 



(30) 



where R{t, s) is given by eq. / f^) . 

The form of the two-point function does not change if we use alt^ instead. 

Ageing systems are often described in terms of Langevin equations of the 
form 

where r]{t) is assumed to be a Gaussian noise with mean zero and variance 
{T){t)r]{t')) = 2T6{t — t'). In the context of Martin-Siggia-Rose theory, it can 
be shown that response functions can be written as a correlator R{t, s) = 
{(j){t)(f){s)) of the order parameter and its associated response field 0. On 
the other hand, to each solution (j){t, r) with positive mass Ai > there 
exists a conjugate solution 0*(t, r) which satisfies the diffusion equation with 
A4 replaced by —Ai and which can be said to have negative mass. 

On the other hand, if we work instead with the wave function ilj{(,t,r) 
as defined in eq. (fT^. we can use the embedding sci)^ C conf^_,_2 the 
form of the two- and three-point functions is fixed by conformal invariance. 
For ageing systems, time-translation invariance does not hold and we must 
restrict to those parabolic subalgebras which do not contain We are 
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interested in the linear response functions 



6h{s) 



(32) 



R3(ti,t2,s) 



6h{s) 



(33) 



Instead, we may also calculate them using the functions tlj{(,t,r) and then 



Proposition 11: |37] For an ageing system with z = 2 and which is invariant 
under agc^, the two- and three-time spatio-temporal response functions may 
be obtained from 



such that the causality conditions ti > s and t2 > s are automatically satis- 
fied. 

This suggests the identification of the response field 



with the conjugate solution of the Schrodinger equation. Explicit expressions 
for these response functions are known. 

Finally, Galilei invariance holds de rigueur only for a vanishing tempera- 
ture T = and in principle all predictions made on ageing systems so far are 
restricted to that case.^ However, constructing the free-field Martin-Siggia- 
Rose action it can be seen through the Wick theorem that the perturbative 
series in T terminates after the first term. 

Proposition 12: j^^j For an ageing Gaussian system with z = 2 and which 
is invariant under agc^ and at a temperature T < T^, the response function 
R{t, s) agrees to all orders in T with eq. / f^) . 

This argument is expected to work in the low-temperature phase, where T 
flows to zero under the action of the renormalization group [2j so that per- 
turbation theory should be applicable. Direct evidence in favour of this is 
given below. The independence of the form of R{t, s) on the (short-ranged) 
initial conditions is all the more remarkable since for instance the correlators 

^At r = 0, the ageing behaviour is determined from the properties of the initial state. 



find 





(34) 
(35) 



(36) 
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do depend on them On the other hand, at T = thermal fluctuations 
and fluctuations of the initial state become simultaneously important. Re- 
markably enough, there exists good numerical evidence that eq. ()29|) even 
holds at T = Tc in the 2D and 3D Glauber- Ising model fl7|. However, a 
two- loop renormalization group calculation of the 0{n) field theory gives a 
numerically small correction to eq. (j29|) [HUj . The understanding of the role 
of temperature on the validity of Galilei invariance in field theory will require 
more work. 

A quantitative test of the prediction is best carried out through the 
consideration of an integrated response function, such as the thermoremanent 
magnetization 



where the system is quenched in a small (random ) magnetic field h which 
is turned off after the waiting time s has elapsed. The magnetization is then 
measured at a later time t. Local scale invariance furnishes, via (j^^ . an 
explicit predictions for the scaling function fM{x)- In certain cases there 
may be sizeable correction terms to this scaling form, see [71 HUl HH E21 EH] 
for details. 

We now consider ageing in the Glauber- Ising model quenched to a 
temperature T < Tc from a fully disordered initial state. Then the exponents 
a = 1/z = 1/2 and Xr ~ 1.26 and 1.6 in 2D and 3D, respectively, are 
known. In figure E^b, the scaling function obtained from large- 

scale simulations, is shown for several values of the waiting time s. In both 
two and three dimensions, a nice scaling behaviour is found and the form of 
the scaling function agrees very well with the prediction from eq. ((221) • 

For z = 2, we have seen in section 3 the importance of Galilei invariance 
as the second building block, besides dynamical scaling, of the local scale 
invariance of systems whose action 5* transforms according to eq. (j23|) . It 
is therefore important to test Galilei invariance directly, which can be done 
by considering the spatio-temporal response R{t, s; r) and testing it for the 
form eq. (jHUj) [52 . We first fix the mass M. by studying how the integrated 
response J^du R(t,u; r) depends on r^/s for a single fixed value of x = t/s. 
Next, a demanding test of the r-dependence of R(t, s; r) can be performed 
by measuring the spatio-temporally integrated response 



MTRM(t,s) = h / duR{t,u) = r^s-'^fMit/s) 



(37) 








(38) 
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Figure 4: Scaling form of the integrated magnetic response in the Glauber- 
Ising model as a function of a; = t/s below criticality. The symbols corre- 
spond to different waiting times s. The integrated autoresponse is shown (a) 
in 2D at T = 1.5 and (b) in 3D at T = 3. An example of the integrated 
spatio-temporal response in 2D at T = 1.5 and with /i = 2 is shown in (c). 
The full curves are obtained from ()29|1 . After |l52j . 

where is a control parameter. We stress that the scaling function p*^^^ does 
not contain any free non- universal parameter at all As an example, 

we compare in figure data from 2D taken with n = 2 with eq. ()29|1 . 
Besides the expected scaling, the functional form of the scaling function 
neatly follows the prediction. We stress that the position, the height and 
the width of the maximum of p*^^-* in figure I3J:; are completely fixed. Similar 
results have been obtained for other values of /z and in 3D as well. This 
provides strong evidence that eq. (jHUj) is exact, at least in this model 
This is the first time that Galilei invariance in an ageing system has been 
directly confirmed. It is remarkable that Galilei invariance, which after all 
was used with a generator which takes the form found for a free particle, is 
confirmed for a theory which is certainly not a free-field theory. 

The prediction ()29|) for the autoresponse has been confirmed in several 
physically distinct exactly solvable systems undergoing ageing, including sev- 
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eral variants of the kinetic spherical model, the voter model and Brownian 
motion, see [U HU 113 EH ES] and references therein. Where applica- 
ble, these models also confirm the form ()30|) of the spatio-temporal response. 

Taken together, these confirmations suggest that should hold inde- 
pendently of 

1. the value of the dynamical exponent z 

2. the spatial dimensionality d> \ 

3. the numbers of components of the order parameter and the global sym- 
metry group 

4. the spatial range of the interactions 

5. the presence of spatially long-range initial correlations 

6. the value of the temperature T 

7. the presence of weak disorder 

This provides strong evidence that local scale invariance is indeed realized as 
a dynamical symmetry in ageing systems. 

7 Conclusion 

Motivated by the observation that dynamical or strongly anisotropic scaling 
occurs in many physically relevant situations, we have made the hypothe- 
sis of dynamical scaling the starting point of our study. Since conformal 
transformations may be seen as a sort of space-dependent local scale trans- 
formations, we have asked ourselves whether an analogous extension might 
be possible for any given value of the dynamical exponent z. Indeed, we 
have shown that infinitesimal transformations of this kind can be explicitly 
constructed. It turns out that there exist two distinct realizations of these 
local scale transformations (Type I and Type II) which lead to very different 
physical consequences. 

These purely formal considerations appear to have a bearing on physics, 
since the predictions of local scale invariance for some two-point functions 
could be explicitly verified in concrete physical models, which describe either 
Lifshitz points (Type I) or ageing systems (Type II). Several of these models 
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do not reduce to a free-field theory. We have also indicated in the text some 
of the open problems which remain and which we hope might be resolved at 
a later time. 

We hope that in the future, the principle of local scale invariance might 
become useful in order to derive hitherto unexplained properties of dynamic 
(or strongly anisotropic) phase transitions. For example, is there a way to 
predict the values of the exponents which we have used here as externally 
given parameters ? 

This work was supported by CINES Montpellier (projet pmn2095) and 
the Bayerisch-Franzosisches Hochschulzentrum (BFHZ). MH thanks the Cen- 
tro de Ffsica Teorica e Computacional (CTFC) of the Universidade de Lisboa 
for warm hospitality, where this work was finished. 
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